arXiv:hep-th/9808155vl 25 Aug 1998 


MS-TPI-98-19 
August 1998 
hep-th/9808155 


Extended Iterative Scheme for QCD: 
the Four-Gluon Vertex 


L. Driesen and M. Stingl 

Institute for Theoretical Physics I, University of Munster 
D-48149 Munster (Westf.), Germany 


Abstract: We study the self-consistency problem of the generalized Feynman rule ( nonperturba- 
tively modified vertex of zeroth perturbative order ) for the 4-gluon vertex function in the framework 
of an extended perturbation scheme accounting for non-analytic coupling dependence through the A 
scale. Tensorial structure is restricted to a minimal dynamically closed basis set. The self-consistency 
conditions are obtained at one loop, in Landau gauge, and at the lowest approximation level (r = 1) 
of interest for QCD. At this level, they are found to be linear in the nonperturbative 4-gluon co¬ 
efficients, but strongly over determined due to the lack of manifest Bose symmetry in the relevant 
Dyson-Schwinger equation. The observed near decoupling from the 2-and-3-point conditions permits 
least-squares quasisolutions for given 2-and-3-point input within an effective one-parameter freedom. 
We present such solutions for Np = 2 massless quarks and for the pure gluon theory, adapted to the 
2-and-3-point coefficients determined previously. 




1 The generalized Feynman rule T^'fy 


The present paper continues, and brings to a provisional stage of completion, the determination of 
generalized Feynman rules in an extended perturbation scheme for QCD (!]], designed to account for the 
strongly nonanalytic coupling dependence of correlation functions through the renormalization-group 
invariant mass scale A. The generalized rules, denoted r[ r ’ 0 l, are proper vertex functions of zeroth 
perturbative order ( no power corrections in the coupling g 2 , as indicated by the index 0 ), but with a 
nonperturbative A dependence which in turn is approximated systematically at a level characterized 
by the index r. In contrast to the ordinary Feynman rules, it is a nontrivial self-consistency problem for 
these extended rules to reproduce themselves in the integral equations for vertex functions. However, 
the divergence-related mechanism operative in the self-reproduction jlj] does ensure that formation 
of Fool’s remains rigorously restricted to the small number of vertices corresponding to the ordinary 
Feynman rules - the superficially divergent vertices. In the companion paper [;2|], this self-consistency 
problem was set up and solved for the vertices with up to three external legs, on the lowest nontrivial 
approximation level (r = 1, and one loop) of interest for QCD. 

Below, we complement this study by examining the highest superficially divergent amplitude, 
the 4-vector (4-gluon) vertex r 4 y [|3j. There are several features that set this vertex function apart 
and motivate its separate treatment. One is its sheer kinematical complication, due partly to the 
large number of 6 Lorentz-scalar momentum variables, but mostly to the exorbitantly lengthy tensor 
decompositions implied by its 4 color and 4 Lorentz indices. In constructing the r = 1 approximant 
(section 2 ), one is thereby forced to adopt a theoretically motivated simplification of the full 
tensor structure. Another special feature is that T 4 y represents the lowest QCD amplitude which in 
the nonperturbative context exhibits the phenomenon of compensating poles ( or, in the terminology 
of [||, of negative shadow poles ), first noted in Abelian models in 0. While these poles can be 
inferred uniquely already from the “lower” DS equations for Tsy and V 2 y by suitable residue- taking 
operations, it is now necessary to demonstrate that on the level of the 4-point equation they are in 
turn self-consistent ( section 3 ). A natural by-product of this demonstration is then a “reduced” 
integral equation for the remainder vertex V 4 y defined by 

lty 0] = -(ci r] ) 2 y,2v -(cP) 2 y,2v - (cf) 2 v, 2 v + V& 0] , (1.1) 

where C{ 2 3 are the shadow-pole terms in the 3 channels of the 4-gluon amplitude. This amounts to 
rearranging the full connected-and- amputated 4-gluon correlation function as 

Tyy = A) + A(> + A 3 + V 4 V, (1-2) 

where A ' 12 3 are “softened”, i.e. one-shadow-irreducible, one- gluon exchange diagrams formed by 
subtracting the shadow-pole term from the ordinary dressed-one-gluon-exchange diagram Ai^ 3 in each 
of the 3 channels. The object V^ v 1 depends on only four Lorentz- scalar variables, and establishing 
the self-consistency conditions for its nonperturbative coefficients (section 4) is simpler than for the 

full r 4 y. 

The latter conditions will exhibit a problem that presumably besets all approximate treatments of 
vertex equations with three or more external lines : if such vertices possess a symmetry among their 
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external legs such as Bose symmetry, their Dyson-Schwinger (DS) equations will in general not display 
this symmetry manifestly, and any approximation to them will produce unsymmetric terms. Forcing 
self-reproduction of a symmetric input by requiring these terms to vanish leads to overdetermination 
of the dynamical conditions. If this does not happen to be counteracted by underdetermination 
tendencies, as was the case for the r = 1 three-point conditions derived in ||], one must settle for a 
“compromise” solution in the sense of least-squares error minimization (section 5).- Technically, the 
self-consistency conditions will be derived at the r = 1 and one-loop (/ = 1) levels, and in Landau gauge. 
The latter again provides considerable simplification because the two ghost vertices at / = 1 remain 
perturbative, and because in this gauge a closed DS system can be written for the amplitudes with 
only transverse (if any) gluon legs, whose tensorial complexity is significantly lower. A disadvantage, 
of course, is that the self-consistency ( or lack of it ) for statements concerning longitudinal-gluon 
amplitudes, such as Slavnov-Taylor identities, cannot be checked directly from such a system. 

The 4-gluon vertex also occupies a special position in that it belongs both to the class of super¬ 
ficially divergent “basic” vertices having their own ( ordinary or extended ) Feynman rules, and to 
the class of amplitudes capable of developing bound-state ( here, glueball ) poles - a phenomenon 
otherwise restricted to the superficially convergent higher vertices. In the present paper we are exclu¬ 
sively concerned with the generalized Feynman rule, a zeroth-order quantity, which like its ordinary 
counterpart contains no information as yet about bound states. The latter, if present, arise from an 
infinite partial resummation of the higher-order, quasi-perturbative corrections 

[g{vfY • rj lv p \{k}, A, V ) (p = 1,2,3,...). (1.3) 

As the notation recalls, these have their own nonperturbative A dependences, which may be com¬ 
plemented by an inversely logarithmic “perturbative” dependence when reparametrizing the quasi- 
perturbative series in terms of a s (k 2 ) rather than g 2 (n). Upon ladder or similar resummation, these 
may build up a pole in some Mandelstam variable ( the squared sum of some subset of the momenta 
{k} ) at some multiple of A 2 representing the square of a bound-state mass. Such a pole is accompa¬ 
nied by a pair of residue functions, each representing an infinite power series in g 2 starting at least 
with p = 1. It is therefore an example of nonperturbative A dependence embedded in an otherwise 
quasi-perturbative term ( starting at least at order p = 2 ) of the amplitude. These characteristics 
should prevent confusion between bound-state poles and the poles which in the present scheme arise 
as rational approximations to branch cuts in the squares of single-external-line momenta, and in ze¬ 
roth quasi-perturbative order. The 4-gluon function is unique in exhibiting both phenomena. It is yet 
another matter ( and not touched upon in this article ) that in a confining theory the calculation 
of S-matrix elements will require a still larger set of generalized Feynman rules, comprising hadron- 
to-quarks and hadron-to-gluons bound-state vertices for the outer corners of S-matrix diagrams - 
vertices of an intrinsically hybrid nature, since their dependence on the squares of single-quark mo¬ 
menta would have to be determined and approximated consistently with the generalized Feynman 
rules for the internal lines and vertices, whereas their dependence on the total virtuality of the bound 
state is governed by “resununed-perturbative” bound-state dynamics. 
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2 Structure of the r = 1 approximant 


Due to the large number of possible tensor structures of a four-gluon amplitude ( e.g., 43 independent 
Lorentz tensors for a fully transverse amplitude, each combined with 8 independent color tensors 
for SU(3)c), one is forced, at the present stage, to adopt some theoretically motivated restriction to 
a smaller subset of tensors. We briefly describe such a restricted form, partly repeating material from 
appendix A of Q to keep the discussion self-contained. The amplitude 

(v [ lv(Pi,P 2 ,P 3 ,P 4 )) ab ^ d (Pi +P2 +P 3 +P 4 = 0) (2.1) 

is expanded over an 18-member set of color <S) Lorentz tensors, of which 15 are linearly independent. 
They are formed from building blocks 

(z = 1...6, J = 1...3), (2.2) 

where the color tensors C^\ fourth-rank objects over the adjoint representation of SU(3)c, are given 
by 


C a bcd ~ 5 c d , 

C 'abed fabnfedm 


z'-k 2 ) _ X X 

u abed ~ u ac<Jdb, u abed 

/^<( 5 ) _ f f /~i(6) 

^abed JcLcnJdbni abed 


$ ad $bc i 
fadnfben •> 


(2.3) 

(2.4) 


in terms of the SU( 3) structure constants f a b c . Only five of these are linearly independent, since 
C'( 4 ) + C^ = 0 by the Jacobi identity. As compared to the most general structure, this 

set omits three tensors containing a symmetric SU(3) structure constant d a b c ■ The Lorentz-tensor 
building blocks L q-) are the three dimensionless ones, 

Lpj u/ = 8 kX 8 ^ u , L k ^ v = 5^5 vX , L^ v = 5 KU 5 XtJj , (2.5) 

which are linearly independent. In obvious ways, all these fall into crossing triplets, where second and 
third members in each triplet arise from the first member by the crossing operations, 


: (2,3,4) - (3,4,2), (2.6) 

= (2,3,4) - (4,2,3), (2.7) 


the s channel being by convention the (1 + 2) «-> (3 + 4) channel. Here 2 stands for the set (j> 2 , A, b ), 
etc. 


The set (2.2-2.5) is distinguished in that it constitutes the smallest dynamically closed subset of 
the full tensor structure, i.e.: 

(i) It closes under Bethe-Salpeter (BS) iteration in each of the three channels of the amplitude, as 
exemplified by the s-channel, color-space multiplication, 


( C® • C (i) ) 


abed 


— ^2 C a bef 

e ,/=1 


0 C U) 


efed 


( 2 . 8 ) 


and closes also under iteration of DS interaction diagrams, provided the 3-gluon vertices involved have 
fabc color structure only. We emphasized earlier ||| that with 4-gluon color dependence restricted 
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to the five-dimensional basis (|2.3|/ |2.4j) , it would be inconsistent to retain d a bc color structure in the 
3-gluon vertex, since the tensors omitted on the 4-gluon level would also contribute to that structure 
via the hierarchical coupling in the 3-gluon DS equation. 

(ii) It closes under the crossing operations connecting the three channels, as given by (^1]) and 

s 


(iii) It contains the perturbative zeroth-order vertex, namely, 

■p(0)perf _ 3 „ia\ / _ N 1 


4 V 


( 


L (2) - L { 3) 


) + 2 ( C& - C{6) ) (~ 2L W + L ( 2 ) + L (3) ) ' 


(2.9) 


( Without this property, there would exist a still smaller dynamically closed set comprising only the 
tensors ([D^), but maintaining the perturbative limit calls for inclusion of ( |2.4| ) ). All these statements 
are easy to verify by direct calculation. 

For the choice ( |2.5| ) of Lorentz tensors, additional motivation comes from the observation that 
these are the only ones accompanied by invariant functions of zero mass dimension. The tensors 
omitted in ( |2.5| ) contain two or four powers of momentum, and are thus associated with invariant 
functions of mass dimensions —2 or —4. Together with the fact that their zeroth-order forms must 
have at least one power of A 2 in each term, this severely restricts possibilities for these functions to 
receive the divergent loop contributions required for self-reproduction, and it appears reasonable then 
to start with a basis omitting such tensors. 

Invariant functions for the V4y amplitude, we recall, have zeroth-order rational structure only with 
respect to the variables Pi 2 ■■■P4 2 , since the entire zeroth-order rational structure with respect to 
s, u, t variables is borne by the compensating-pole terms of eq. (1.1). Since our earlier demonstration 
of this applies, strictly speaking, only to the color-octet channels, we check in appendix A of this 
paper that in fact it holds generally. Thus at order r = 1, all invariant functions involve the same, 
fully symmetric denominator 

A [1 ’°] = ( pi 2 + u' 2 A 2 ) (p2 2 + 'u'jA 2 ) (p 3 2 + '4'A 2 ) (p 4 2 + ?4'A 2 ) . (2.10) 

Numerator polynomials in A 2 and the squared momenta, of mass dimension 8 , must conform to the 
restrictions of naive asymptotic freedom. As seen in || in the example of the 3-gluon vertex, additional 
restrictions arise from the requirement that “softened” ( one-shadow-irreducible ) one-gluon exchange 
diagrams, which represent the physical one-gluon exchange mechanism, again should not generate 
higher-than-perturbative degrees of divergence when inserted in loops ( for an example, see the right- 
hand portion of diagram (A 3)4 in Fig. 10 ): this excludes terms in the rational approximant in which 
any one squared momentum has a net positive power. Thus numerator polynomials should not contain 
s, t, u either, nor should they have terms of type A 2 p\ or A 4 p 4 still allowed by the primary restrictions. 
At level r = 1, all zeroth-order rational structure can then be expressed in terms of the single-pole 
quantities 

A 2 

IT = —-ttttt ( i = 1,2,3,4), (2.11) 


Pi 2 + u'2 A 2 

and the general r = 1 approximant, with full Bose symmetry among the four legs, turns out to involve 
seventeen dimensionless, real coefficients, each preceding a product of a tensor of type (222) and a 
combination of one or more quantities 
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In the following we use two different representations for . The first has building blocks with 
manifest Bose symmetry but involving linearly dependent color tensors internally; it uses and for our 
purpose defines the 17 independent coefficients. It is thus suitable for the presentation of the vertex 
as a generalized Feynman rule. This representation uses a color tensor 


C ibld = tT (tatdtbtc) + tr {t a t c t b td) 


( 2 . 12 ) 


and its crossing partners C^ U \C^ obtained via ( 2 . 6 ) and (2/7), where t n are the generator matrices 
of SU(3)c normalized by tr (t m t n ) = \5 mn . The latter trace is then used to supply another tensor 
triplet, proportional to ([h3|), 

l(1) (2.13) 


1 


tv(t a t b )ti(t c t d ) = -C^, 


plus crossing partners. SU (3) representation algebra [j| 

C (s) = — [ c (1) + c (2) + c (3) 

12 


gives the relation of ( |2.12| ) to (2.3/2~I|) 

+ - [ C (5) - C (6) 

6 


as 


(2.14) 


with analogous relations for C^ U \C^ T \ In this basis, which has been used e.g. in work on the 
operator-product expansion Q, the vertex with four transverse gluon legs is written 


■■■Pi) 


pCTTUJ 

- abed 
17 


= t p ^ Pl )t^(p 2 )t^( P3 )t^(p,) 


{(r'r‘)2; + 

k 71=1 


kX/iv 'l 
abed 




(2.15) 


the t’s now denoting transverse projectors in the external momenta. Each of the building blocks W n 
comprises a color tensor ( |2.12| ) or ( 2.13 ) plus crossing partners, a combination of Lorentz tensors (p.5|) , 
and a combination of rational-function elements ( 2.11 ). Each therefore has at least one overall power 
of A 2 and vanishes in the perturbative limit, A —> 0. A full listing of these building blocks is relegated 
to appendix B.l; here we write out just one of them for illustration: 


»TV...j<4 2 ;A 2 ) 


c (s) ( 

+ CIO ( 

+ C< T > ( 


^(3) - 2T(i) + L(2) 
L (l) ~ 2L (2) + £(3) 

L(2 ) - 21/(3) + ^(1) 


II1II4 + 112113 
II1II2 + II3II4 
nin 3 + 1I4II2 


2 (nin 2 +113114) + nin 3 +114112 
2 (nin 3 +114112) + iii 1I4 + n 2 n 3 
2 (nin 4 + IT2II3) + ihiE + 1I3II4 


(2.16) 


Note that these 17 building blocks are linearly independent of each other; the linear dependence 
of the basis tensors C^ S,U,T ^ affects only their internal structure, which upon rewriting in terms of 
independent tensors would lose its manifest crossing ( and therefore Bose ) symmetry. However, as a 
representation of the total V 4 y, ( 2.15 ) is based on dependent color tensors. 

The second representation is more suitable for the self-consistency calculation of the zeroth-order 
vertex. To permit tensor-by-tensor comparison of the two sides of a DS equation, such a representation 
must employ linearly independent color tensors, and consequently forgo manifest Bose symmetry. 
However, since the way a 4V vertex appears in a DS interaction diagram singles out one of the three 
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channels ( compare, e.g., the 4V vertex in diagrams (A 3)4 or ( 6 * 2)4 of Fig. 10 ), it is useful to maintain 
partial symmetry properties within one channel. A basis of this kind, adapted to the s channel, has 
color tensors 


c*( A ) 


C {B) 

= C (2) 

+ c* (3) , 

= 

C (5)_ C (6), 

(2.17) 

C (C) 

= c (2) - 

c (3) , 

cW 

= c (4 ) 1 

l=-c^ 

- cw) , 

(2.18) 

and Lorentz tensors ( 

in D = 4 - 

- 2e Euclidean dimensions ), 





L (0) = ^ 

L (i)’ 

L (+) 

2 

= 2 ^ 

1 

rZ 

+ 


(2.19) 




L (-) 

= 2 ^ 

- £(3)) • 


(2.20) 


The latter have been chosen to display, under s-channel tensor multiplication, the projector-like prop¬ 
erties, 

L t7 L ?4" = i-C (m, rc — +,0, —). (2.21) 

If we denote by the symbol 


{ tti 2 , 7r 34 ; ir s } ( each tt = +lor — 1 ) 


( 2 . 22 ) 


the class of four-point amplitudes having parities 7 ri 2 , ^ 34 , 7 r s under the 1 2, 3 <-> 4, and (1,2) <-» 

(3,4) permutations respectively, then of the fifteen product tensors formed from ( |2.17|/ 2.18| ) 

and ( |2.1£ / 2.2C ) there are eight in class {+, +; +} and seven in class {—, —; +}. In the decomposition 

(F^ 0] ( P1 ...P4))27 = ^(Pl)^(P2)t T,i (P3)t W,, (P4) 


' abed 

v V' r>(») P [i,o], 2 2 . A 2 i 

X Z^ ^ abcd L (j) *i,j \P 1 ■•■P4 >A ), 

i=A...E 0) — 


(2.23) 


the invariant functions F t j therefore come in two types, 


F, 


[1,0] 


(M')e{+,+;+} 


F, 


[1,0] 




= Vi,j, 0 + Vi,j, 1 (III + II 2 + n 3 + II4) + j 2 (III + II2) (II3 + II4) 

+ 4 (1I1II2 +113114) + r/j.ys [nin 2 (n 3 +114) -1- (iii +112)113114] 

+ Vi,j ,6 (IIin2ll3ll4) ; 

= Vi,j ,3 (Hi — H2) (1I3 —114). 


The perturbative limit (2.9) fixes the coefficients 77 *°f ( 2.24 ) as 

Vi,j,0 = T &i,E [ (T) 1 )^ 0 ] • 


(2.24) 

(2.25) 


(2.26) 


This representation then starts with a total of 8 x 5 + 7 x 1 = 47 real, dimensionless coefficients. 
With no further restrictions on the latter, it would be adequate for quantities having partial Bose 
symmetry with respect to one distinguished channel, such as a two-particle irreducible kernel. For the 
full V^y amplitude, imposing complete Bose symmetry produces 30 relations between the 77 ’s so they 
can ultimately be expressed linearly in terms of the 17 independent £’s of ( 2.15 ). The listing of these 
linear conversion equations is relegated to appendix B.2. 
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The r.h.s. of a DS equation such as Fig.2 or Fig. 10 below, which singles out the leftmost external 

( although it 


leg in an unsymmetric way, does not even exhibit the reduced symmetry of (2.24/2.2 
does have a permutation symmetry between the three rightmost legs ). Upon using ( |2, 24 / ^25 ) as 
input, the output therefore has still lower symmetry, and enforcing the self-reproduction of a symmetric 
vertex turns out to imply, in addition to 47 self-consistency conditions proper, 7 more conditions for 
the vanishing of terms excluded by the partial Bose symmetry, leading, after conversion of rj 1 s to £’s, to 
a total of 54 linear conditions on the 1 7 independent coefficients. Moreover, one finds that unwanted 
contributions to the eight constants ( [2. 26 ) arise which represent low-r approximation errors in the 
perturbative renormalization constant, a phenomenon already encountered in || for the three-point 
vertices. 

Finally note that with respect to any one of the four pf. the reduced vertex permits a partial- 
fraction (p.f.) decomposition, such as 

r's? 1 = <I(ft 2 ,ft 2 ,ft 2 ) + t ( P 2 + A „ 2 2 „A, ) (ft 2 ,ft 2 ,ft 2 ), (2,27) 

where by omitting pf dependence in the Eq term we have encoded the restriction of no net positive 
powers of pi . Such a decomposition is useful technically in connection with the treatment of the 
compensating poles. 


3 Self-consistency of compensating poles 

It was demonstrated in [j2j that from the lower ( i.e., T 3 y and T 2 V ) equations alone one may already 
infer the existence, and determine the residues, of certain pole terms in T^y with respect to the three 
Mandelstam variables, and that these turn out ( as did similar poles in longitudinal channels in the 
work of refs. j|] ) to be compensating terms cancelling unphysical singularities (“shadow” poles ) in 
the one-gluon reducible terms of T±y. The compensating poles are structures of tree topology which, 
since their internal shadow lines are not propagators of any of the elementary QCD fields, nevertheless 
can be present in the 1PI amplitude. This result is recalled in eq. (0) and in Fig. 1(a), where each 
double-wiggly line, at level [r, 0], stands for a set of r shadow poles in the corresponding Mandelstam 
variable. 

To understand the nature of these structures more fully, we now need to discuss how they reproduce 
self-consistently in the DS equation for T^y itself. As with all zeroth-order nonperturbative terms, this 
will essentially occur through the hierarchical DS coupling, which transfers the building materials for 
the compensating terms down from the higher vertices in the equation. The argument is lengthy, but 
conveys an idea of how, by analogy, such terms establish themselves in still higher vertex functions. 
We therefore present it also in lieu of a general N-legs-amplitudes proof, whose length would be out 
of proportion to the insight gained. ( The reader not interested in the details of this self-consistency 
result, and willing to accept it on faith, may proceed directly to the reduced DS equation for V 4 y 
given diagrammatically in Fig. 10, which is the starting point for sect.4.) 

(1 ) The equation for T^y is written diagrammatically in Fig.2, in a “hybrid” form in which 
interaction terms on the r.h.s. have not been resolved down to the level of proper vertices: the T' 
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amplitudes appearing there are amputated functions 1PI in the channels defined by their right-hand 
external legs, but otherwise still contain 1PR terms. 

To keep technical complications to a minimum, we invoke all available simplifications: we disregard 
fermions, i.e. the term (E )4 of Fig.2, and consider zeroth-order self-consistency only at one loop, where 
term (D) 4 does not yet contribute, and in Landau gauge, so that the ghost term (B) 4 may also be 
omitted. We are then dealing, first, with terms (C)4, (C 1 ) 4, {C") 4, involving the four-gluon function 
itself in the form of three T' 4V amplitudes, each one-gluon irreducible in the channel defined by its 
two right-hand external legs. The structure of such a T' 4V as following from Fig.1(a) is recalled in 
Fig. 1(b): it has one compensating ( negative shadow ) pole in the distinguished channel, and “softened” 
( one-shadow irreducible ) one-gluon exchanges in the two other channels. 

( 2 ) Second, we have term (A )4 involving a five-gluon amplitude T(y, which by definition is one- 
gluon irreducible in the four 2V <-> 3V channels accessible through its three right-hand external legs, 
namely the “horizontal” channel 

(5, 6 ) <- (2,3,4) (i = 1), (3.1) 

and the three “tilted” channels 

(2,3) (4, 5, 6 ) (t = 2), (3,4) ^ (5, 6 , 2) (i = 3), (4, 2) <- (3, 5, 6 ) (i = 4). (3.2) 

The integers i assigned refer, by convention, to a numbering i = 1... 10 of the ten 2V 3V channels 
of a five-point amplitude. Thus T^ v has those one-gluon reducible terms removed that would produce 
a one-gluon reducible graph when introduced into term (A) 4 . A representation of T§ v suitable for our 
purpose is given diagrammatically in Fig.3; it is again hybrid in that it involves T' 4V amplitudes in 
addition to fully 1-gluon irreducible functions T^y and Tsy. The derivation of this representation is a 
technical matter that we relegate to appendix C. ( Its apparent asymmetry with respect to the (5,6) 
pair of legs is resolved by realizing that there exists an equivalent representation with the roles of T^y 
and T' w interchanged ). 

Insertion of Fig.3 into the (A )4 term of Fig.2 then decomposes that term into a piece (A) 4 .p 
involving the 1-gluon irreducible T$y, plus six pieces with triangle-graph topology. We do not draw 
this decomposition separately. 

( 3 ) Now take residues in the T4y equation of Fig.2 with respect to the squared momentum, p \ 2 , of 
the leftmost external leg. This is entirely analogous to the residue comparison performed in [|| for the 
T 3 y equation, and leads to the analogous conclusion: the zeroth-order pole terms in p\ 2 of the r.h.s. 
can arise only from the term (A^p, and the Tsy in that term must therefore contain zeroth-order 
poles (si + « 2 nA 2 ) _1 in the Mandelstam variable of channel ( |3.1| ), 

«i = ( -P5-P6 ) 2 = (P 2 +P 3 +P 4) 2 , (3.3) 

which equals p\ 2 by momentum conservation. Note that terms with an inverse polynomial in p\ 2 
produced by the loop integration of term (A) 4 .p, if any, would by the standard formulas of dimensional 
integration have to arise from convergent parts of the integral, and therefore would not trigger the self- 
consistency mechanism for zeroth-order terms; they would be part of the first-order , quasi-perturbative 
correction - a remark indeed applying to all instances of generation of zeroth-order terms. 


The residues of Tsy at those poles must be proportional to the residue functions £+ ( P 2 ,P 3 ,P 4 ) (1 < 
n < r) in the p.f. decomposition of analogous to eq. ( [2.27 ). On the other side they must be 
proportional to the residue functions B^i (/+,/+) (1 < n < r) of the three- point vertex Tsy; this 
follows e.g. from the way the same T$y amplitude enters in the two-DS-loops term of the three- point 
equation ( term (D )3 in Fig.l of || ). For these pole terms of Tsy, the loop of (Aj^r then turns into 
loops of self-energy type already encountered in the self-consistency problem of , which fixes the 
proportionality factors: the internal lines in the zeroth-order pole terms of T$y with respect to si are 
associated with factors 

1 t(pi) 


- ++) = - 


(n = 1... r ) 


(3.4) 


Ur t 2n+1 Pi 2 + Ur,2n A 2 
and therefore constitute minus a shadow line as defined in Fig. 3(b) of Q. 

( 4 ) To demonstrate self-consistency of the three negative-shadow terms in eq.(Ll), we now feed 
the decomposition of Fig. 1(a) into the various contributing terms on the r.h.s. of Fig.2 and verify 
that, upon appeal to information from the “lower” DS equations, those terms cooperate so that their 
sum splits off explicitly the same triplet of negative-shadow terms; after cancelling them from both 
sides of the equation, the remainder then constitutes a DS equation for the reduced amplitude V 4 y. 

The first place to feed in Fig.1(a) is the T^y poles just identified, with their £$ residue functions. 
This turns each of them into a pole with the corresponding reduced E n instead, minus a triplet of 
terms with two shadow poles each, involving in addition to the B\ r ) (n > 1) a set of functions Bm, n 
defined by further decomposition of these, 

A 2 

- 1 h . 

m,n\ 


B [ £(p,q) = < 0 (<?) + £ 


B [r] (q) 

™ n V'i / 


(3.5) 


\p 2 +U rt 2n^ 2 / 

Invoking now the full Bose symmetry of Tsy, we see that this vertex has a shadow-poles structure 
given, in a condensed notation, by 


F5V = Fj;y 
10 

- E 


(3.6) 


2=1 

+ E{E s t!sB B anSL r|s t'b’ 

k =1 m,n 

where i enumerates the 10 channels, while k enumerates the 15 different 2+1+2 partitions of the 5 legs. 
This structure is shown diagrammatically in Fig.4. Note that the symmetrization is achieved with 
only 15, not 30, of the latter terms: since each of them has shadow lines in two different channels, these 
are enough to supply the triplet of contributions mentioned before to pole structure in all of the 10 


Mandelstam variables. The amplitude V$y defined by eq. (|3.6|) will be referred to as the reduced 5-gluon 
vertex. As a superficially convergent vertex with all tree-like structures removed by construction, it 
consists only of superficially convergent loops. Since these, as emphasized repeatedly, do not support 
the self-reproduction of zeroth-order terms, it is an exact statement that 


y[ r fi\ _ n 

v 5V ~ U 


(all r' 


(3.7) 
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( 5 ) The next place to feed in the assertion of Fig.l(a), or its consequence of Fig.1(b), is the six 
terms of triangle topology, arising from (A) 4 of Fig. 2 upon insertion of the 1-gluon-reducible terms 
in the 2nd and 3rd lines of Fig.3. Together with eq.( |3.6| ) and after some regrouping of terms, this is 
equivalent to using a new representation of the T~ v amplitude given diagrammatically in Fig.5. We 
again relegate details of its derivation to appendix C and note only that this representation, apart 
from the fully one-shadow-irreducible gluon-exchange graphs in its 2nd and 3rd lines, has a shadow-line 
content isolated in four terms: the negative-shadow term in the first line, which is identical to the i = 1 
term in the 2nd line of eq. (|3.6| ) and will play a special role by supplying in Fig.2 the singularities of 
V 4 v in its Pi variable, and the 3 negative terms in the last two lines, which involve a 3-gluon-1-shadow 
auxiliary amplitude E n as defined in Fig. 7. 

( 6 ) Lastly, we should introduce Fig.1(b) into the (C) 4, (C 1 ) 4, and (C") 4 terms of Fig.2. This will 
turn each of them into a graph having a Viy instead of Tiy, plus a pair of one-shadow-irreducible 
gluon-exchange terms giving equal contributions and for which the factor of -7 therefore gets cancelled, 
plus a negative term with one shadow line, as displayed in Fig. 8 . ( Note that an ordinary gluon line 
connecting to at least one bare vertex, such as the Y^y ert in Fig. 8 , has no shadow content ). 

( 7 ) In the crucial step of the argument, we now refer back to the three-point vertex equation, 
maintaining for consistency the simplifications analogous to those of step (1) above, and again feed 
the assertion of Fig.l into the T^ v term of that equation ( term (^ 4)3 in Fig.l of [2|], where the T^ v is 
denoted T' s ). We then take residues with respect to the squared momentum of one of the two external 
legs other than the “distinguished” leg entering from the left ( one of the variables pf or p| in Fig.l 
of Q ). Note that this kind of residue information stood unused so far, since in Q we only exploited 
residue comparison with respect to the distinguished leg. The result is the diagrammatic relation of 
Fig.9, a DS equation of sorts for the partial amplitudes B n (n > 1) of Tsy. ( If we do not invoke the 
simplifications, this relation of course gets additional ghost-loop, fermion-loop, and two-gluon-loops 
terms with their own E-type amplitudes, as alluded to in the last line of Fig.9 ). 

But this relation is precisely what is needed to turn the 3 negative shadow-line terms which arose 
in step ( 6 ) from the (C) 4 through [C") 4 terms of Fig.2, plus the 3 negative E terms that arose 
from term (A )4 upon insertion of the E graphs of Fig.5, into the 3 negative one-shadow terms on the 
r.h.s. of Fig. 1(a). We have therefore attained our goal: we have shown that upon introducing the 
decomposition of Fig.1(a) with its 3 compensating-pole terms, either directly or through its immediate 
consequence of Fig. 1(b), into wherever a four-point building block occurs in the interaction terms of 
Fig.2, these interaction terms cooperate so that their sum splits off explicitly those 3 compensating 
terms again. Note that in so doing we appealed to residue information not only from the lower ( 
3-point and 2-point ) DS equations, but now also from the 4-point equation itself, which led us to 
infer the shadow-reducibility structure ( |3.6D of the next higher, 5-gluon vertex. 

We may now cancel the 3 compensating terms on both sides of the T 4 y equation and thus establish 
the DS equation of Fig.10 for the reduced vertex Fiy, which is the main result of this section. 

(8) As a by-product, we may rewrite eq. ( |3.6| ) for the shadow content of Tsy. To each of 
the 10 terms with negative signs in the second line, we may add a triplet of terms of structure 
X). m T 3 yDB m S m B m , with one ordinary-gluon and one shadow line, such that the resulting quartet of 
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terms form precisely the set of compensating poles needed to cancel the four shadow lines present in 
the one-gluon-exchange diagram {Tsy DT^v}i of the full connected-and-amputated 5-gluon amplitude 
T$y ( one along the D propagator line, the other three hidden, by Fig.l(a), in the T±y vertex ). The 
same 30 terms, now regrouped into 15 suitable pairs, may then be subtracted from the 15 terms in 
the third line of ( |3.6| ), and it is straightforward to check with the aid of Fig. 6 (b) that each resulting 
triplet then forms the set of compensating poles needed to cancel the one two-shadows term and two 
single-shadow terms present in a double-gluon-exchange diagram {TsyDTsyDr^y}^ of T$y. Since the 
total operation has not changed the equation, we conclude, with a view to eq. ( |C.1| ) of the appendix, 
that 


r 5V = V^y + [ complete set 

of compensating poles for all 25 one-gluon-reducible terms of T§y 


(3.8) 


Thus the final result for the next higher amplitude is again simple: the full T$y is completely one- 

shadow irreducible, consisting as it does of a V$y with full extended irreducibility, plus 25 “softened” 

( i.e. one-shadow irreducible ) but still one-gluon reducible diagrams. 

( 9 ) Finally note that term (Al )4 of the reduced V 4 y equation of Fig. 10, upon taking residues 

at its horizontal-channel shadow poles and invoking once again the F^? self-consistency conditions, 

r r i _ .___ 

instantly reproduces the Ej lt terms with n > 1 of the p.f. decomposition (2.27). The remaining, 

further-reduced equation for the partial amplitude of (2.27) at r = 1 is what we shall work with 

in practice; it should be kept in mind that this equation still has a term comprising the regular-at-poles 

remainders of (Al )4 of Fig. 10. 


4 Self-consistency conditions 

4.1 The overdetermined system 

A glance at Fig.10 shows that the reduced DS equation for V 4 T, or for its i?o,T component, is linear 
in the four-point unknowns, though still nonlinear in the “lower”, 2-and-3-point vertices. This sim¬ 
plification results from (i) the restriction to a one-loop calculation and (ii) our use of the ordinary 
DS system without Bethe-Salpeter resummation. The self-consistency equations for the ( coefficients 
of eq. ( [2.15 ) will thus form a linear system but with a matrix ( determined by the divergent parts 
of diagrams (Al) 4 , (A 3 ) 4 , and (C2)± of Fig. 10 ) and with inhomogeneous terms ( determined by the 
divergent parts of the other diagrams of Fig. 10 ) depending nonlinearly on the 2-point and 3-point 
coefficients. 

As discussed in Q, the zeroth-order self-consistency system as a whole exhibits a scaling property , 
arising from the scheme- insensitive character of the self-reproduction mechanism. It implies that the 
system only determines ratios of the nonperturbative coefficients to certain fixed powers of one of 
them. The latter, by convention, was chosen to be x±, a coefficient of the 3-transverse-gluons vertex. 


We need to maintain this convention for the 4-gluon quantities; thus for the parametrization (2.15) 
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we introduce rescaled parameters 

A 2 = MA 2 , Ci = Ci/\xi\ ni (i = 1...17), (4.1) 

where the integers rii are, 


ni 

= n 7 

= n i3 = 

= 1, 

n 2 

= n 3 

= 74 6 = 

?4 8 = 749 = 7412 = 74 i4 = T4l 7 = 2 

744 

= n w 

= 7415 

= 3, 

74 5 

= nu 

= 7416 

= 4. 


(4.2) 


All RG-invariant masses are obtained in terms of A 2 rather than A 2 . 

For the actual evaluation of the zeroth-order terms of diagrams (Al) 4 , (A 3 ) 4 , (C 2)4 of Fig.10, 
one uses parametrization ( |2.23 ) with its independent tensors and 77 coefficients. In contrast to what 
happens ( at level r = 1 ) in the 3-point systems, the immediate self- reproduction of the entire 
p | 2 -singular portion E\ F of eq. ( |2.27 ), achieved by extracting the pole term at p 2 = —U 2 A 2 = —U 2 A 2 
from diagrams (Al )4 and appealing to the 2 -gluon self-consistency conditions, is found not to provide 
substantial relief for the overdetermination caused by the lack of manifest Bose symmetry. After 
reimposing the symmetry by expressing all ry’s in terms of the ( coefficients of the fully Bose- symmetric 
form d), one is still faced with 54 linear equations for these 17 basic coefficients, with matrix 
elements and inhomogeneities depending nonlinearly ( though polynomially ) on the lower-vertex 
coefficients. The complete system is listed in Appendix D. To convey an impression of its structure, 
we write here just one typical equation: self-reproduction of the coefficient f?A,o ,4 hi the invariant 
function ( 2.24|) associated with the tensor structure gives, after conversion to C’s by means 

of appendix B.2 and multiplication by /?o and u 2 , the condition 

2 V“3 ^3J Cl + 2 i 8^ 3 ~ U3X4 / ^ + 9 ( ( 8 ~ W u 3 “ 9^3X4 I C 3 
+^(u 3 a:i)C4 + j(^u 3 2 + u 3 x^j ( 6 + | + Po) u 3 2 ~ 


15 

- — x 3 - N f z 3 


/" 

x (C8 — 2(9 + 2£l2 — Cl4 + C17) + U3 ( 

9 (- 7 }U 3 2 x 3 2 + 4 u 3 xix 3 x 4 + (^u 3 - x 2 )x 2 - 2xix^x^j 


(C10 - C15) 


+ N f 


u £\ 2 2 

w 3 ) 3 


Z A {w 3 ZlZ 3 - W 3 Z 4 + Z!Z 5 - Z2Z4) • 


(4.3) 


Here Nc = 3 has been used. The Np factors identify contributions from the quark-loop diagrams 
(£71)4, ( j E 71 / )4 of Fig.10. Note that for the quark-vertex coefficients Zi, the simplified numbering of 
eqs. (4.6) of [§] has been employed. 

As far as we are aware, no presently existing mathematical software tool for nonlinear algebraic sys¬ 
tems is capable of dealing directly with a system of this size and complicated, partially overdetermined 
structure. The only strategy currently practical for dealing with the total system is to take advantage 
of the “near decoupling” of the 4-point system: of the 17 ( coefficients, only two combinations, Z\ and 
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Z‘i , enter the 2-and-3- point system (through some of the 3-gluon vertex conditions). One then treats 
these as additional unknowns in the lower-vertices system and, after invoking the scaling property, 
determines them to within an one-parameter freedom. Subsequently, their definitions are appended 
to the linear system for the 4-point £’s as two extra conditions, bringing the total to 56 equations: 
this represents only a minor increase in the anyway massive overdetermination of that system. The 
extra equations are, in rescaled form, 

15(30-4) - Z, - |) (4.4) 

^(3C2 + 3C3-C8-C 9 ) = ^2 = ^§, (4.5) 

oZj X 

where the Z\ : Z 2 on the r.h.s. stand for the values produced by the lower-vertices system. The 
structure of the four-point system then is, 

17 

J2 M ik({u},{w}-,{x},{z})(k = bi({u},{w}-,{x},{z}) (4.6) 

k= 1 

( i = 1... 56 ), 


where the 56 x 17 matrix M, as well as the 56-component inhomogeneity b, depend nonlinearly on the 
rescaled coefficient sets {h}, {if:} of the gluon and quark two-point functions and on the sets {x}, {z} 
of the gluon and quark three-point functions. Exceptions are the 55th and 56th rows of M, given by 
the l.h. sides of eqs. (L4 / 4. 5|) , which contain only pure numbers. 

For a quasi-solution in the least-squares sense, one minimizes the quadratic deviation between both 
sides of (IO). The conditions 


r) 56 

E ( E - bif } = 0 (n = 1...17) 


(4.7) 


i=l 


lead, in standard fashion, to a linear system with the 17 x 17 matrix M T M, 


17 


E(M T M) nfc a = ( M T b\ 

k=1 


( n = 1... 17 ). 


(4.8) 


Although the structure of M is too complicated to be handled analytically, we have checked numerically 
that in the physically acceptable parameter range where all 2-and-3-point coefficients are real, the 
matrix M T M is always invertible. ( The calculations have been performed using the MAPLE V 
computer-algebra system ). As a measure of the deviation one may consider 


X = 


\ 


1 

56 


56 

E (J2 M ^k-bi 


(4.9) 


i=1 


as compared to a typical r.h.s. or l.h.s. of eqs. (L6). Also, the quantities ^ 1,^2 when recalculated 


from the quasi-solution according to eqs. (L4/L5) will be different from their input values from 
the lower-vertices system, and the differences may serve as rough indicators of the overall degree of 
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difficulty in satisfying zeroth-order self-consistency requirements from ( ordinary ) DS equations with 
the rather simple structure of the r = 1 set of approximants. If they are large, one may alternatively 
try to enforce conditions (4.4/[Q|) exactly by assigning a large weight w (e.g. w = 1000) to the i = 55 
and i = 56 terms of eq. (T7). The Z\ 2 from the two systems will then match, but the quality of the 
54-term remainder solution will in general deteriorate. 


4.2 Typical solution for N F = 2 

The entire system of course shares the effective one-parameter freedom from the “decoupled” lower- 
vertices system. In the presence of fermions, this was parametrized [j||] by the coefficient w\. We only 
consider here the physically most interesting parameter range where not only all vertex coefficients are 
real, but also the values of the self-energy coefficients u, w imply the presence of a complex-conjugate 
pole pair in both the gluon and the fermion propagator. Over this range, as noted earlier, most of the 
other coefficients of the system do not vary substantially. 

Table 4.1 gives least-squares values of the rescaled coefficients ( corresponding to the “typical” 
solution of table 5.1 of Q, which at | w\ | = | W 2 |= 0.6749 ( the signs of w\,w 2 do not affect the 4- 
point system ) is in about the middle of this parameter range. The two solutions shown are obtained 
either without extra weighting {w = 1) or with exact enforcement {w = 1000) of the 55th and 56th 
conditions. The resulting coefficients are generally rather large in absolute value, given the fact that 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Cn (W = 10 3 ) 

56.22 

646.67 

-513.11 

611.96 

962.88 

271.12 

-168.19 

-803.84 

871.27 

Cn (W = 1) 

-181.48 

-98.04 

-113.83 

297.85 

-2747.9 

-37.55 

172.27 

-319.20 

35.75 

11 

10 

11 

12 

13 

14 

15 

16 

17 


Cn (W = 10 3 ) 

64.77 

438.55 

-202.04 

-510.04 

84.67 

-561.86 

-7481.9 

209.80 


Cn (W = 1) 

1318.4 

-1757.5 

-26.14 

76.28 

190.68 

367.79 

5446.0 

-317.7 



Table 4.1. Coefficients C n for Np = 2 at typical value w\ = 0.67 


from their very context one would expect them to be numbers of order unity ( give or take an order 
of magnitude ). This may partly be attributed to the choice of basis - the factors of |, |^ in the 
tensors (p,13 / pT4 ) employed in the W n building blocks - and to the rescaling (4.1) : it is conspicuous 
that the large coefficients mostly have large exponents n* in ( |4.2| ). Although these effects may account 
for one to two orders of magnitude, several coefficients still remain implausibly large. We interpret this 
as the direct result of the considerable “pressure” generated by the overdetermination of the system. 

The typical 2-and-3-point solution of [p] implies values of 


Zi = 157.9, Z 2 = 156.2 ( N f = 2) 


(4.10) 


for the parameter combinations of eqs. (44/44^) respectively. When not enforced ( i.e. at w = 1), 
their values recalculated from the C' s of table 1, —335.95 and —165.07, are strongly mismatched and 
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even of the opposite sign. It is then hardly surprising to find that this solution has x ~ 317 in the 
sense of ( |4.9| ) - of the same magnitude as, or even slightly larger than, a typical \bi\ in (^1]), which 
is exemplified by ( 4. 10| ) and of the order of a hundred. Enforcing the values ( 4.10| ) would then be 
expected to strongly alter the solution, which is indeed what happens - several of the coefficients even 
change signs, and y deteriorates further to ~ 349. This is in marked contrast to the purely gluonic 
case discussed in the next subsection, and implies that in the presence of ( massless ) fermions the 
4-point solution at the r = 1 level is probably a poor one. The presence of the massless-quark DS 
loops evidently makes the total system at low r much harder to satisfy; in particular these loops seem 
to work in the direction of forcing a mismatch at the (2+3)-point-to-4-point interface that can be 
alleviated only gradually at higher r. 


4.3 Typical solution for pure-gluon theory 

For pure Yang-Mills theory (Np = 0) the typical solution, now parametrized by the gluonic coefficient 
X 3 as discussed in j2j] and taken at X 3 = 1, leads to the 4-point <C’s given in table 4.2. 


11 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Cn (W = 10 3 ) 

0.814 

6.439 

-4.733 

-57.35 

350.05 

-3.057 

8.672 

-30.03 

0.931 

Cn (W = 1 ) 

-0.432 

9.285 

-4.789 

-54.834 

255.08 

-2.512 

6.398 

-21.50 

1.275 

n 

10 

11 

12 

13 

14 

15 

16 

17 


Cn (W = 10 3 ) 

96.02 

-328.14 

-0.060 

3.354 

-4.212 

13.74 

-121.67 

-9.060 


g 

II 

1 - 1 

56.75 

-110.07 

-0.518 

2.595 

0.461 

-19.15 

124.0 

-8.744 



Table 4.2. Coefficients Cn for pure-gluon theory at typical value X 3 = 1.0 


Two features deserve comment. First, the coefficients overall are of a distinctly more plausible order 
of magnitude, entirely understandable from the above-discussed simple mechanisms; in particular, the 
three coefficients of sizes more than a hundred are precisely those with = 4 in (4.2). Second, and 
more remarkably, the input Z 12 values from the lower-point system, 


Zi = -2.92, Z 2 = 16.04 (AY = 0) 


(4.11) 


are matched reasonably well by the four-point quasisolution even for w = 1 , i.e. without being enforced 
exactly: when recalculated they come out as -3.60 and +15.80 respectively. Moreover, the x value of 
~ 0.11 is now only about ten percent of the typical | 6 j| of (|4.6D, which in this case is of order unity. 


Exact enforcement of the values ( 4.11 ) then leads to no appreciable deterioration in y. This gives one 
confidence that for the pure Yang-Mills system the least-squares, four-point quasisolution does work 
reasonably well even on the r = 1 level of approximation - better, in fact, than one would expect in 
view of the still very simple and rigid structure of approximants at this level. 
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5 Conclusion 


We have verified that self-consistent determination of a generalized Feynman rule for the highest 
superficially divergent and kinematically most complex QCD vertex is possible in principle, provided 
the overdetermination resulting from lack of manifest Bose symmetry in the relevant DS equation is 
dealt with by a least-squares procedure. Within the well-defined and clearly visible limitations of our 
calculation, available indicators suggest that the least-squares solution at level r = 1 is good for the 
pure gluon theory but quantitatively poor in the presence of massless quarks. 

Among the limitations, the one most obviously in need of improvement may be the use of an 
“ordinary” DS system: starting from BS-resummed dynamical equations for N > 3 would presumably 
shift more important physical effects into the lower loop orders, and would provide partial ( though 
not complete ) relief from the pressures of overdetermination. Similar effects can probably be expected 
from use of a tensor basis larger than our “dynamically minimal” one; this would moreover allow d a b c 
color dependence in the three-gluon vertex to be treated consistently. Going beyond the r = 1 and 
l = 1 levels of approximation would have the algebraic complication rising steeply, as usual in iterative 
schemes for QFT, but is presumably the only consistent way of better satisfying perturbative limits 
and other desirable secondary conditions. The above results do, however, give one confidence that the 
method as a whole is viable, and produces enough nontrivial features to render these improvements 
worth pursuing. 
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Appendix 


A Absence of non-compensating zeroth-order poles 


• f y* 0l 

Here we ask whether the four-gluon generalized Feynman rule F^, may contain, in addition to the 
’’compensating” poles in its three crossed color-octet channels as identified in sect. 2.2 and made 
explicit in eq. (1.1), still other poles of zeroth perturbative order in its Mandelstam variables. Such 
poles would be absent in the perturbative limit A —> 0, so their residue factors would be proportional 
to at least one power of A. In e.g. the s channel, such a pole would be of the form 

A^> T (pi,p 2 )A^(p3,P4) , A -n 

P 2 + 6A 2 [ } 

with a dimensionless residue function (or vector of functions) carrying the color and Lorentz quantum 
numbers of the channel considered. This contrasts with the case of the compensating pole, where 
residue comparison in the T 3 y equation forces the residue factors to be proportional to B n 's - quantities 
of mass dimension +1 which by themselves need not vanish as A —> 0 since in their defining equation 
( eq. (2.14) of ref. (^] ) they are already accompanied by a A 2 factor. 

The Bethe-Salpeter amplitudes A<h must satisfy a well-known normalization condition [^], which 
in a condensed notation, and again in the most simplified form omitting ghost and fermion terms, 
reads 


+ 1 90J 


4/ M, (fl?d DD> ) A ‘ 1, 

11 M>{DD)(fL K ,)(DD)M> 


p2=-bA 2 


= -2 P p . 


pz=-bA 2 


(A.2) 


This is depicted in Fig. 11. Fulfillment of this condition by purely zeroth-order quantities is possible 
only if the loop integrals on the l.h.s. diverge and thereby trigger the divergence-related 1 /g 2 mecha¬ 
nism. With the $ functions behaving at worst like constants at large loop momenta q, and with the 
easily checked behavior of 


d 


■( DD) = 0(q ~ 5 ) 


d 


■K s = Oiq- 1 ) 


(A.3) 


d pu ' ' ' ’ dP^ 

(note that terms in K s behaving like constants at large q' s, such as the Y^ pert term, drop out after 
differentiation), it is however clear that the integrals are convergent, and do not provide the j factors 
necessary for the mechanism to work. Thus poles of type ( |A-l| ) in r]y°^ are ruled out. 

Note that the ’’compensating” poles , by contrast, can satisfy ( ]A.2| ) as zeroth-order quantities, 
since there A<h is replaced by a B n ( compare eq. (2.21) of j|] ) which at large q behaves like q l , and 
therefore provides divergences on the l.h.s. of ( A.2 ). 

Note also that if one does not insist on the being quantities of zeroth perturbative order, ( |A.2| ) 
may be satisfied, at suitable eigenvalues — bh? of P 2 , by quantities which have no overall A factor 
but are power series starting at least at first order in g 2 . This, as emphasized in the introduction, is 
the case for ordinary bound-state poles, which arise from a Bethe-Salpeter partial resummation of the 
perturbation series for Y^y- 
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B Parametrizations of the reduced four-gluon vertex 


B.l Building blocks for the symmetric parametrization 

To list the building blocks w\\ ’°^ of the fully Bose-symmetric V 4 y representation ( 2.15 ), we use the 
following combinations of single-pole quantities ( 2 . 11 ): 


G(i) ip 1 ■ ■ ■ Pi ) — III + II2 + II3 + II4 ; (B.l) 

G( 2 s){Pi ■ ■ -Pa) = H1II2 + 113114 , G(2u)(Pi ■ ■ -pI) = H1II3 + Ibn^ , 

g( 2t){pl ■ ■ ■ p!) = nin 4 + n 2 n 3; (B.2) 

G( 3 ){pI ■■■Pi) = nin 2 n 3 + 111112114 + nin3ii4 + n 2 ri3n4; (B.3) 

G(4)(.Pi • • • Pi) = IhrblTiILt . (B.4) 


Here G*(i), G( 3 ), and G( 4 ) are Bose-symmetric in themselves, while G 2s , G 2u , and G 2 t form a crossing 
triplet. In terms of these, 


<’"1 : 

— | C (S ' ) [1/(3) — 21/(4) + 1/(2)] 



+( perms, s —> u, s —> t ) j Gp) ; 

(B.5) 

wf °1 : 

= { C^ s ) [I/ (3) — 2L (1) + Z/ (2) ] 



+( perms, s —> u, s —> t ) j (G( 2s ) + G^ u ) + G^t)) ; 

(B.6) 

wf' 01 : 

= | [1/(3) - 21/(4) - 1 - 1/(2)] 



X [G( 2 t) - 2 ^( 25 ) + G {2u )} } + ( perms, s -> u, s -> t ) 5 

(B.7) 

W M . 

= { G (s) [L( 3) - 2L (1) + L( 2) ] 



+( perms, s —> u, s —> t ) | G( 3 ) ; 

(B.8) 

, 

= | [1/(3) — 21/(4) + t'(2)] 



+( perms, s u, s —> t ) | G( 4 ) ; 

(B.9) 

W M . 

= [^(l) + L (2) + ^(3)] 



X [G( 2 t) - 2G( 2s ) + G (2u )] } + ( perms, s -> u, s ->• t ) ; 

(B.10) 

W™ : 

= {-C (1) [1/(3) - 21/(4) + £(2)] 



+( perms, s —> it, s —> t ) | Gp) ; 

(B.ll) 

<’°1 : 

= {~C (1) [1/(3) - 2L (1) + 1/(2)] 



+( perms, s —► u, s —► t ) | (G( 2s ) + G( 2u ) + G( 2 t)) ; 

(B.12) 

: 

= [L(3)-2L(4)+L(2)] 



X [G( 2 t) - 2G( 2s ) + G( 2u )] } + ( perms, s -> u, s -> t ) ; 

(B.13) 
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wM . 
*' 10 

- {-C (1) [L (3) - 2L (1 ) + L (2) ] 



+( perms, s —>u, s —> t ) | ; 

(B.14) 

W/M : 

= [-^(3) “ 2L (1) + -^(2)] 



+( perms, s —> u, s —> t ) | G(4) ; 

(B.15) 

W/M : 

= {^C' (1) [L(i) + L(2) + L( 3) ] 



X [G( 2 t) ~ 2G (2s ) + G'( 2u )] } + ( perms, s -► u, s -> t ) ; 

(B.16) 

<’ 01 : 

= ^P (1) + C (2) + C (3) ] [L (1) + L (2) + L (3) ] G (1) ; 

(B.17) 

wt" 1 1 

= ^P (1) + + C®] [L (1) + L (2) + L (3) ] 



X [^(2t) + + ^(2u)] j 

(B.18) 

< 01 : 

= ^p (1) + + tfP)] [L (1) + L (2) + L (3) ] G (3) ; 

(B.19) 

< 01 : 

= ^P (1) + + C^] [L (1) + L (2) + L (3) ] G (4) ; 

(B.20) 

W^f ] 1 

= ^P (1) + ^ (2) + C (3) ] { [^(3) - 2T (1) + L (2) ] G {2s) 



+( perms, s —> u, s — > t ) [ . 

(B.21) 


B.2 Relation between parametrizations 

To obtain the expressions for the 7 ] parameters of ( 2.24| ) and ( 2.25 ) in terms of the 17 £ coefficients, 
one rewrites the Wn'^ building blocks ( as listed above ) in eq. ( 2.15| ) in terms of the linearly indepen¬ 
dent color tensors ( j2.17 / |2T8| ) and the s-channel-adapted Lorentz tensors ( 2.19 / T2C| ) and compares 
coefficients with eq. (2.23). The ensuing relations are, 


VA, 0,1 = 

3 A 3 A 

" 2 C? + 2 Cl3 

(B.22) 

VA, 0,2 = 

3 a 3 a 3 a 3 a 3 a 3 a 

-4C 3 - 2^ 8 ~ 2 < ’ 9 + 2^ 12 + 2^ 14 + T^ 17 

(B.23) 

VA, 0,4 = 

2C3 - 2^8 + 3< ^ 9 — 3 Ci2 + 2^ 14 _ 9^ 17 

(B.24) 

VA, 0,5 = 

3 A 3 A 

- 2 < ’ 10 + 9 815 

(B.25) 

VA,0,6 = 

3 A 3 A 

" 2°i + 2 Cl6 

(B.26) 

Va,+, 1 

1 u 1 . 

- 2C7 + 2C13 

(B.27) 

^?A,+,2 = 

1 A 1 . 1 A 1 A 1. 1 

4C3 + 2 ’8 + 2 ’9 + 2^ 12 + 2^ 14 ~ 4^ 17 

(B.28) 

^,+,4 = 

“ 2 + 2 ~ ~ ^ 12 + 2 ^ 14 + 2^ 17 

(B.29) 

Va,+, 5 — 

2 C10 + ^ Cl5 

(B.30) 
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VA,+, 6 

= 

1 . 1 . 

2 ^ n + 2 ’ 16 





(B.31) 

VA-, 3 

= 

“ I Cl7 





(B.32) 

Vb, o,i 

= 

3 A 3 A 

4 C7 + 2 Cl3 





(B.33) 

VB, 0,2 

= 

3 A 3 

- + 4 C 8 - 


- 

3 

+ 

3 . 3 

2 C 14 + 4 C 17 

(B.34) 

VB, 0,4 

= 

3 a 3 > 3 

2 Cs + 4 Cs + 4 

C 9 + 

3 

2 

A 3 

Cl2 + 2 

3 

C14 - ^Cl7 

(B.35) 

0,5 

= 

3 l. 3 

^ClO + ^Cl5 





(B.36) 

VB, 0,6 

= 

3 a 3 

4 C 1 I + 7 J C16 





(B.37) 

? 7 B ,+,1 

= 

“ 4 C7 + 2 Cl3 





(B.38) 

VB,+,2 

= 

J~Y 

co 

1 

h-> 

</~> 

00 

+ 

00 | 1— 1 

C 9 “ 

1 

4 

Cl2 + \ 

C14 - 4 C 17 

(B.39) 

VB,+, 4 

= 

-^Ca - 4 Cs ^ 

4 C 9 + 

^C !2 + 

^C!4 + ^Cl7 

(B.40) 

T?B,+,5 

= 

1 1 

- 4 ClO + ^Cl5 





(B.41) 

? 7B,+,6 

= 

“ 4 Cll + q C16 





(B.42) 

VB,-,3 

= 

3 9 3 

4 C 3 + gCo - 4 

C17 




(B.43) 

VC, 0,3 

= 

- ^Co - |Cl 2 





(B.44) 

VC,+,3 

= 

^C 0 - ^12 





(B.45) 

Vc,-, i 

= 






(B.46) 

VC,-,2 

= 

3 . 3 

“4 Cs + 8 C9 





(B.47) 

VC,-A 

= 

“ I Cs “ I C9 





(B.48) 

VC,-,5 

= 

- l^o 





(B.49) 

VC,-,6 

= 

“I Cn 





(B.50) 

VD, 0,3 

= 

9 9 

4 Cs + 2 Ce 





(B.51) 

VD,+,3 

= 

3 A 3 

~ 4 C 3 + 2 < ’ 6 





(B.52) 

VD,-, 1 

= 

i<* 





(B.53) 

VD,-,2 

= 

^C 2 - |Ca 





(B.54) 

VD,—,4 

= 

3 A 3 u 

2 C 2 + 2 C 3 





(B.55) 
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VD,-,5 = 

5 <* 

(B.56) 

VD,-, 6 = 


(B.57) 

Ve, 0,1 = 

-§* 

(B.58) 

VE,0,2 = 

3 A 3 > 3 A 

- 2 82 _ 4 83 + ^ 86 

(B.59) 

VE,0,4 = 

3 3 

- 2 + 2 _ 

(B.60) 

VE, 0,5 = 

-!«* 

(B.61) 

VE, 0,6 = 

~l Cb 

(B.62) 

VE,+, 1 = 


(B.63) 

VE,+, 2 = 

1. 1 A 1 A 

2C2 + 4 83 + 

(B.64) 

VE,+,A = 

2C2 - -(3 - C 6 

(B.65) 

VE,+, 5 = 

1 „ 

2 

(B.66) 

VE,+, 6 = 

2 Cs 

(B.67) 

VE,-, 3 = 


(B.68) 


C Representations of the amplitude T§ v 

Here we sketch derivations for two different representations of the partially irreducible 5-gluon ampli¬ 
tude Tgy needed in the argument of sect. 3. 

( 1 ) To establish the representation given in Fig.3, start from the full connected and amputated 5- 
gluon amplitude ( off-shell 5-gluon T matrix ), T§y. It has 25 one-gluon-reducible terms ( dressed tree 
diagrams ): one of structure T3y DF^y, reducible along a single internal gluon line, for each of the 10 
different 2+3 partitions ( or 2V 3V channels ) of the 5 legs, and one of structure T^y DF^y DT^y, 
reducible along two different internal lines, for each of the 15 different 2+1+2 partitions of the 5 legs. 
In the condensed notation adopted for eq. (|3.6| ), we have 


Tsv = F 5 y 
10 

+ (C.l) 

2—1 
15 

+ VI {F3y DF 3 yDF 3 y } fc . 

k= 1 


Of these one-gluon-reducible terms, 13 are reducible in at least one of the four channels (3.1/3.2), and 
therefore by definition are to be excluded from the Tf ]V amplitude, leaving the latter to consist of the 
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fully one-gluon-irreducible piece ( proper vertex ) Tsy, plus 12 dressed-tree diagrams still reducible in 
other channels - 6 from the second and 6 from the third line of (C.l). 


The 6 one-gluon exchange terms from the second line of (C.l) are those in which one of the two 
legs no. 5 and 6 entering from the left in Fig.3 connects to the T3y, while the other connects to the 
T^y. The 3 terms in which leg no. 6 connects to a T3y constitute the second line of Fig.3. The 3 
terms in which leg no. 5 connects to the T3y may each be taken together with a suitable pair from 


among the 6 terms retained of the third line of (C.l), in such a way as to replace their F^’s with 
T^y's whose primes refer to the channels defined by their external-line pairs. ( Remember that, by 
definition, T^ v = T^y plus 2 one-gluon exchange terms in the channels other than the distinguished 
channel to which the prime refers ). This establishes the third line of Fig.3. 

Of course, in this regrouping the roles of legs 5 and 6 may be interchanged, leading to a repre¬ 
sentation equivalent to Fig.3 in which the T^/s and T^’s have their roles exchanged. This explains 
the apparent asymmetry of Fig.3. ( 2 ) The representation of T^ v given in Fig.5 is established by 


( i ) applying the decomposition of eq. ( |3.6| ) to the T^y in the first line of Fig.3, ( ii ) inserting the 
decompositions of Fig. 1(a) and Fig. 1(b) into the Tyy and the T' AV terms of Fig.3 respectively, and ( in) 
performing decompositions of the remaining ordinary gluon lines into “softened” one-gluon exchanges 
and shadow-pole terms as indicated in Fig.6. ( As for Fig.6(b), it lumps together for simplicity two 
decompositions, corresponding to its first and second lines, that actually occur separately ). The six 
last terms of Fig.3 then give rise to: 

(a) 6 softened one-gluon-exchange graphs of the type of term (B) of Fig.6, shown in the second 
line of Fig. 5, 

(b) 6 graphs with two softened one-gluon exchanges each, of the type of term (E) of Fig.6, shown 
in the third line of Fig. 5, 

(c) 6 one-shadow-line graphs as in term (C) of Fig.6, which cancel 6 of the 10 negative terms 
coming from the second line of the Tsy decomposition ( |3.6|) , leaving just those four of the 10 negative 
terms that possess a shadow line in one of the four channels (|3.1|/ [3.2D , 

(d) 12 graphs with two shadow lines each of the type of term (H) of Fig.6 - 9 from the T^y and 
3 from the T' 4V terms - with minus signs, which cancel 12 of the 15 terms from the third line of eq. 
(3.6), leaving just those three of the 15 terms that have one shadow in the “horizontal” channel (p.l|) 


and the other in one of the three “tilted” channels (3.2), 

(e) 18 graphs with one softened-gluon exchange and one shadow line each: 12 with minus signs 
from both the Tyy and T' 4V terms, and 6 with plus signs from the T^ v terms alone, and which actually 
cancel 6 of the minus-sign graphs, leaving a net count of 6 terms with one softened-gluon and one 
shadow line and with minus signs. 

Now the 6 terms surviving from (e), the 3 two-shadows terms surviving from (d), and 3 of the 4 one- 
shadow terms from (c) may be lumped together in three groups by defining the auxiliary amplitudes 
E n as in Fig.7; they then give the three graphs in the last two lines of Fig.5. The one remaining one- 
shadow graph with minus sign from (c) appears explicitly in the first line of Fig.5, which is therefore 
fully established. 
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D Listing of V^y self-consistency equations 

In the following, the 54 self-consistency conditions ( |4.6| ), multiplied by /3qu 3 to rid them of ubiquitous 
factors from the self-consistency mechanism, are listed according to the C® ® L^) tensor combination 
[i = A... E, j = 0, +, —) in whose invariant function they arise. For the quark-vertex coefficients Z{, 
the simplified numbering of eqs. (4.6) of [^] is employed. 

D.l Ca~Lo equations: 


-U3X3 ^5 £2 + C3 + 8^6 — 6^8 + 12 C9 — 12Ci2 + 6^14 — 6^17^ — — 773X4(4 
+ (J/3 oM3 + (C10 - C 15 ) + (^U 3 X 1 - y#3l 3 - N F U 3 Z 3 ) (Cll ~ Cl6 ) 

= — 27^3X3 + 54U3X3X4 — 27X4X5 


(D.l) 


7^3X3 Cl + - 7^3X4) C2 - (^Poul - ^ul + ^ 773X4) C3 + ^3X4(4 

+ \~^ u 3 + ^3X4) C 6 + (2^° W 3 -Y ~ ^“3X4) (^(s - 2( 9 + 2(i2 - C14 + C17) 

f99 15 Ar \/ \ 

+ M3X1- — 773X3 — N f U 3 Z 3 J (^ C10 — Cl 5 J 

27 2 2 9 2 2 

=-— 77 3 X 3 + 36773X4X3X4 + -773X4 — 9X2X4 — I8X1X4X5 

, r / 773 \ 2 / 2 2 2 2 2 2 \ 

+ Np ( - - 777304 ^ 3^4 - - W 3 Z 4 + - 04 - 74-35 - - 2 2 2 4 

V 7773 / V 3 3 3 3/ 


(D.2) 




tol CO 


/ \ / 3 9 \ 9 

ulx 3 [ 2 Cl - 3 C 7 + 3 Ci3 J - 9 U3X4 (2 + Po + - U3X4 J C3 - 2 U3X1 C4 

+ (y0 U 3 ~^ U 3 x ijC6 + (^Aj«3 + — U 3 X 4 J (2Cs + 2Cg - 2Cl2 —2^14 - C 17 ) 

+ {^U 3 X 1 - y«3^3 - N F U3Z^j (ClO - C 15 ) 

45 2 2 2 2 

=-— U3X3 36113X1X3X4 + 27443X4 — 9x’ 2X 4 — 18x1X4X5 

Ar /U 3\ 2 /1 1 9 1 1 o\ 

+ N F ^ — J y ~ - W3ZxZ3Z4 + - W3Z4 - - Z X Z4Z 5 + - 2 2 2 4 ) 

(D.3) 



(D.4) 
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\ U &3 ( 4 Cl - C 7 - Cl 3 ) + ( | u\ - 2 U3X4 ) C2 - ( J Po ul + | 113 + ^ U3X4) Cs 
- ^ U3X1 C4 + ul + y w 3 x 4 ) Ce - Ad u\ - i u§) ( 2 Cs + 2 Cg 

+ 2(12 + 2Cl4 — C 17 ) — (—^ 3^1 - -U 3 X 3 - - NFU 3 Z 3 ^ ^ClO + C 15 ) 

=-— U3X3 + 12 U3X1X3X4 + -1(3X4 — 3X2^4 — 6X1X4X5 

, , / U3 \ 2 / 1 1 o 1 1 o\ 

+ 7 V F - - 1 X 3 2 i 2 3 2 4 - - 10324 + - 212425 - 7T 222 4 

V vj-i / V 9 9 9 9 / 


(i u 3 “ ^“3X4) Cl - -J U3X1 (l 4 C 2 - 13 Ce) 

+ (21(3X1 - ^ U3X3 - i N f u 3 z 3 0 (3 - (^Ad« 3 “ 7^113X4) (Cr + C13) 
— (— W3X1 — -U 3 X 3 — -N F U 3 Z 3 ^ (Cs - 2 Cg — 2 C 12 + Cl4 + C 17 ) 


= — — U3X3 + 6^32:1X3 + 5W3X1X4 — 3 xfx 5 — 6 xiX 2#4 

, r /w 3 \ 2 / 2 2 2 2 2 2 \ 

+ M —] 9^2124 + - 1 X 32123 - - 212224 + - 2i2 5 J 


^t«3 - \ u 3X 4 ) Cl - £ U 3 X 1 ( 16 C 2 - 11 Ce) 

(y« 3 Xl - “7 K 3 X 3 - ^N F U 3 Z 3 ^( 3 - - jju§) (C 7 + C 13 ) 

(-g-W — —“ 3 X 3 — -N f u 3 Z 3 ^ (2Cs + 2Cg + 2 C 12 + 2Ci4 — Ci7) 


= — - M3X3 + 61(3X1X3 + 81(3X1X4 — 3X1X5 — 6X1X2X4 

/ K3 \ 2 / 1 1 2 1 1 2 \ 

+ —J (9^32124 - -1x32123 + -212224 - -24 2 5 J 



















ool CO 


ulx 3 ( 8 Ci - 2 ( 7 + Ci 3 j - 9 u 3 x 4 (2 + ( ^ Po uj -— u 3 x 4 ) C 3 - ^ u 3 X! C4 

/27 2 63 \ /3 2 9 \ /3 2 9 \ 

" V16 U 3 + ~8 U3X4 J* 5 \4 “ 3 “ 4 M3:r4 J + [g P° “ 3 + 4 M3:r4 J & 

/27 15 1 r \ /3 2 27 2 9 \ A 

- (^yW 3 £i - Yq U3X3 - 2 NpU3Z3 ) ^ 10 + V 4^° m 3 “ 16 ~ g 113X4 ) 8 12 

/3 . 2 27 2 9 \ /99 15 Ar \ 

(y Po U 3 ~ ^ u 3 + g w 3^4J Cl4 - [-M 3 I 1 - y U 3 x 3 - N F U 3 Z 3 J Cl5 

/3 2 9 \ * 

(y Po W 3 + g W 3 Z 4 J Cir 


45 2 2 27 2 2 

=-— it 3 2 3 + 36 1 / 30 : 10 : 3 : 1:4 + — u 3 x 4 — 9x 2 x 4 — I8xix 4 x 3 

( u 3 \ 2 / 1 1 9 1 1 9 \ 

+ N f J ^ - w 3 z 4 z 3 z 4 + - w 3 z 4 - - zxz 4 z 5 + - z 2 z 4 J 


27 ,9 / , „ \ /45 15 Ar \ „ 

— 1 / 30:4 Ci - ^ 213X4 ^ 2 C 2 ^ Ce J - (^— 22 - 3^1 -— ^ 32:3 - N F u 3 z 3 J ( 3 

Q /fl 1 1C 1 x / x 


jPo uf 


Mi 


15 

u 3 x i - 77 ? 


y U 3 X 3 - \ N F U 3 z 3) ( Cs + Cq ) 

/27 15 Ar \ / A A \ /3 . 2 27 2 \ 

(^yM32;i - Y U3X 3 ~ N FU 3 Z 3 J ^ (,12 + Cl4 J - (y U 3 ~ \Q U 3 ) 

15 Ar \ . 

— u 3 x 3 - Np u 3 z 3 J C 17 


C 13 


/ 15 

^9«3^i- 5- 


= — - u 3 x 3 + 18u 3 xfx 3 + 27 u 3 x\X 4 — 9x 4 x 3 — 18x±x 2 X4 


/ 1/3 \ 2 / 2 2 2 2 2 2 \ 

+ Np ^ — j (y 1 C 3 2 i 2 4 - - W 3 z 4 z 3 + - 210224 - g 2325 J 




(27 2 9 \ „ 9 / A \ / 9 15 1 T \ 

V16 “ 3 - - U 3 X 1 ^ 8 C 2 - "C 6 j + 1*3*1 - J^ U 3 X 3 - 2 NfU 3 Z 3 J <3 

(3 . 2 9 n /27 15 1 \ 

- ^Po «3 - -U3X4JC7 - U 3 Xi - y^W 3 X 3 - -^N F U 3 Z 3 j Cs 

/81 15 1 Ar \ „ (9 15 1 r \ 

+ VIG^ 3 ^ 1 ~ 32 11,3X3 ~ 2 NfU3Z3 J ^ + V 2 U3Xl ~ ~16 113X3 ~ 7 ^ n fu 3 z 3 J C 12 

/3 o 27 2 9 \ (99 15 Ar \ 

- [^P0 U 3 _ ^2 M 3 + g “32:4J C13 - - y U3X3 ~ n fu 3 z 3 J Cl 4 

(27 15 1 r \ , 

- ^—M3a;i - Y^ U 3^3 - — Np u 3 z 3 J C17 

= — ^7*3X3 + I8U3X1X3 + I8W3X1X4 — 9 x ^ x 3 — 18x1X2X4 

at ( u 3 \ 2 ( 1 1 2 1 1 O \ 

+ IVf J ^ - g w 3 ziz 4 + - w 3 z 1 z 3 - - z 4 z 2 z 4 + - z 4 z 5 J 

(D.17) 


D .5 Cb~L + equations: 



(D.18) 
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J U 3 X 3 ( 14 Cl - 14 C 7 + 13 Cl 3 ) - (^ 7«3 + 7^32:4) C2 + (7; A) u\ + ~~ U 3 X 4 ^ £3 
\ U3X1 C 4 + u\ + ^ 1132:4) Ce + Po u\ + 7 - 1x32:4) ( Cs + Cq ) 

+ (y 1X32:4 - Y^ u 3 X 3 - ^ N F U 3 Z 3 ) (ClO - 2(45) 

- (^/^0 U 3 + — 1 X 32 : 4 ) (Cl 2 + Cl4^ — (^2@° U 3 + ~^ u 3 x a) Cl7 

17 2 2 11 o o 2 

=- 2 ~ U 3 X 3 + 12 1X32:42:32:4 + — U 3 X 4 — 6X2X 4 — 62:42:42:5 

.. /«3\ 2 / 2 2 2 2 2 2 \ 

+ ^ { — ) { - 9 + g *> 3*4 - g - 1^425 + g 2 2 2 4 J 

(D. 19 ) 

-^1X32:3(16(4 - 14 C7 + 13 ( 13 ) + (J«3 ~ 2 IX324) (,2 - ( j^Poul + i 1(32:4) C3 

7 . /15 2 11 \ (1 a 2 3 2 7 \ . 

- - 1x32:4 C4 + u 3 + Y 1X3X4 ) ^ + (4 ^0 Ms “ 8 Us + 4 U3:r4 / 

/I . 2 3 2 7 \ A /29 5 1 Ar \ A 

- [^PoU 3 - -u 3 - -U3;r4jC9 + — U3X1 - — 1X32:3 - - -/Vf W 3 2 3 J ClO 

/l ,, 2 15 2 13 \ „ /I 2 15 2 13 \ . 

+ [^PoU 3 - —u 3 - y 1X32:4^(12 ~ P0 U 3 - — U 3 + — 1X32:4J Cl4 

/43 5 1 Ar \ A /I 2 3 2 7 \ v 

- V"gT M3:ri _ 8 U3X3 ^ 3 NfU 3 Z 3J Cl 5 + y^P 0 U 3 - g U3 - — U 3 X 4 J (,17 

11 ^ 

xx OO u o 9 

= —— 1x32:3 + 12 M32:i2;32;4 + — u 3 x 4 — 3x2X 4 — 6 xix 4 x 3 

„ /w 3 \ 2 /l 1 2 1 1 o\ 

+ 7VF I^J (g 1^32123^4 - g»32 4 + g 2l2 4 2 5 - g 2 2 2 4 J 

(D. 20 ) 


3 







(D. 21 ) 



(D. 22 ) 
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D .6 Cb~L _ equations: 


- ulx 3 ( 4 Cl - 4 Ct + 5 Cl3 ) - ( M 3 + 2 U3X4 ) ( ) + 4 u 3 Cs 

+ | M3X1 ( 4 C 4 - 4 Cio + 5 Ci5 ) - (yy u 3 - y 22324) ( Ce - C12 ) + (Ad'< 

3 \. / 9 2 15 \ /3 2 9 2 3 \ 

+ - u 3 24J Cq - U 3 + y w 3 24J Ci4 - (yA d u 3 - —u 3 + - u 3 x 4 J C 17 

= - 3 M 2 a ; 3 + 314^4 + N F (-^-) ( - ^w 3 z 3 z 3 z 4 + ^ it 

1 1 2 \ 

- - .2124.25 + g 2 2 2 4 J 


(^ #3 + 1^32:4) (C1-C7) + g2i 3 £i (4C2 - 5 Ce - 4 Cs + 5 Ci 2 + 5 C14) 
/45 15 1 r \ /159 45 3 Ar \. 

ly U3Xl ~ 16 113X3 ~ 2 NpU3Z3 ) < * 3 ~ \~T6 U3Xl ~ 32 U3X3 ~ 2 NfU3Z3 )^ 9 

, 9 2 15 \ /51 15 1 r \ . 

V39 “ 3 + y M 324j Ci 3 + [y 113X1 - Yg U3X3 - 2 NfU3Z3 )w 

= 3u 3 XiX 4 + Np ( —) (\ W 3 ZiZ4 - l; W 3 zfz 3 

V vj 3 / V 3 3 

1 1 2 A 

+ - 212224 - - Z 1 Z 5 J 


II 


ool co 








D.8 Cc~L + equations: 


U3X4) (C3 - C17) - ^3^3 ( 14 C 7 - 13 <13) + (J/?0«3 + \ U3X4) Cs 


iWo — — Ur 


+ ^3*4) C9 - (l 4 Cio - I3C15) - Po 1 


13 \ , /15 2 13 \ l. 

+ Y 113X4 ) l ’ 12 _ { 32 m 3 + Y U3X4 J < ’ 14 


0 










■U3®1 (7 Cl + 14 Cs - 13 Cl4 - 7 C 17 ) + (J M 3 + J U 3 X 4 ) C 7 - (y^ “a 2 -! - p M3X3 

,, \ . /29 15 1 . \ . /15 2 13 \ . 

■ IVf U323J Cg + y— U3X1 - — U 3 x 3 - - Np U 3 Z 3 J C 12 - U 3 + — M 3 a; 4 J Ci3 


= 0 



1 











^w 3 £i( 2C3 - 5Ci2 - 5Ci4 - 2 C 17 ) + (^A )«3 - ^ ^ 3 ) C 7 + (yU 3 a;i - ^ u 3 o; 3 
- ^NpUzZzj ((8 + C9) + Y^UgCi 3 = 0 


3 / \ /3 9 3 \ 

- « 3 *l(^ 2 C 3 - 5 Cl 2 - 5 Ci 4 - 2C17J + (^A) U 3 - y ^ u 3 + 2 U3X4: ) ^ 

(57 15 1 r \ . ( 21 15 1 r \ 

+ y—u 3 xi - — m 3 x 3 - -N F U 3 Z 3 ) Cs - \ 16 U3Xl - 32 U3X3 ~ j n f u 3Z 3 J Cq 

/ 9 2 !5 \ t. 

+ l 32 U3 ^ y U3;C4 J Ci3 = 0 


D .10 Cd~Lq equations: 

jj «§Z3 ( Ci + C7 + Ci3 ) + u\ - u 3 x^j C 2 + Po u\ + ^ u% + y u 3 x^ Cs 

9 / \ /9 27 45 \ 

+ g u 3 X\ ^(4 + C10 + C15J + ^o «3 — 32 U3 ~2 U3X4 ) 

(27 2 9 w \ (27 2 27 \ 

" (32 U3 + 8 113X4 ) O + Cl4 J - l^ u 3 + y U3X4 ) ^ 

/27 27 \ 9 

- ( 7^7 u 3 + -J- U 3 X 4) C 12 - g u 3 a:4 C 17 

9 o 2 27 2 « 2 « 

= — — u 3 x 3 + 9u 3 Xia; 3 a:4 — — it 3 x 4 — 9 x 2 X 4 + 9x10:40:5 

+ N f ( - W 3 Z\Z 3 Z 4 + w 3 z\ - z 3 z 4 z 3 + 22-24) 








(t^ m 3 ~ g 113*4) Cl - g« 3 ®i (C2 - Cs + 4 (g + 7 C12 - C14 - C17) 

/27 45 3 Ar \ „ /45 45 

- 113*1 - — 1x3x3 - - N f U3Z3J C3 - ^—113*1-^-113*3 - 3 / 

- (i"» + l U3X< ) - (I“» + s” 314 ) <13 


Nf 11323) C6 


= 91x3X4X3 — — W3X4X4 + 9X1X5 — 9X1X2X4 

+ N f (^) (11132124 - W32423 + 212224 - 2^5) 


D.ll Cd—L + equations: 


i ix§x 3 ( 7 Cl + 7 C7 + 13 C13 ) - (p «3 + g 113*4) C2 - (| A) W3 + m H3 

- y 1x3X4) C3 + ^113x1 (7C4 + 7 C10 + 13 Ci 5 ) + (|/?o ii3 - p 113 + 1 1x3x4) Ce 
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Figure 1: (a) Decomposition of 4-gluon vertex into reduced vertex and negative-shadow (“compen¬ 

sating” ) poles, (b) Decomposition of partially irreducible 4-gluon amplitude T' into reduced vertex, 
compensating pole, and softened 1-gluon exchanges. 
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Figure 3: Representation of partially 1-gluon-irreducible 5-gluon T matrix T^ v . T§v denotes the 
proper 5-gluon vertex. An equivalent representation has the roles of T^y and T' iV interchanged. 
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Figure 4: Decomposition of proper 5-gluon vertex into reduced ( 1-gluon- and 1-shadow-irreducible 
) vertex V^y, negative one-shadow terms in the 10 channels, and positive two-shadows terms for the 
15 (2+1+2) partitions of the external legs. 
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Figure 5: Another representation of the amplitude. The E amplitudes in the last two lines are 
defined by Fig.7. 
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Figure 6: Definition of (a) single and (b) double 1-shadow-irreducible (“softened”) gluon exchanges, 
marked by dots, through extraction of shadow terms from ordinary single and double gluon exchanges. 
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Figure 7: Definition of 3-gluon-1-shadow auxiliary amplitudes S. 



Figure 8: Decomposition of the (C )4 term of Fig.2, obtained by use of Fig. 
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+ { ghost, fermion, and 2-loop terms } 


Figure 9: DS-like integral equation for partial amplitudes B n (n > 1) of 3-gluon vertex T^y. The 
auxiliary amplitude E in the second line is defined in Fig.7. 
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— { ghost terms not contributing to Oth order at £ = 0, / = 1 } 
+ { 2-loop terms not contributing to Oth order at l = 1 } 

+ { V .5 terms not contributing to Oth order at any r.l } 


Figure 10: DS equation for reduced 4-gluon vertex V^y, obtained from Fig. 2 upon using Figs. 5, 8, 
and 9. 
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Figure 11: Diagrammatic representation of Bethe-Salpeter normalization condition. The partial 
differentiations with respect to P ^ act on the portions in dashed brackets. 






